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CONTROLLED  DISPERSION  OF  THE  EXTENSIONAL  WAVE  IN  COMPOSITE 
ISOTROPIC  AND  TRANSVERSELY  ISOTROPIC  RODS 

INTRODUCTION  AND  BACKGROUND 

OBJECTIVE  AND  OVERVIEW 

The  objective  of  this  work  is  to  investigate  the  dispersion  of  extensional  waves  in 
composite  isotropic  and  transversely  isotropic  rods.  In  particular,  we  are  interested  in  the  material 
and  physical  parameters  that  control  the  maximum  and  minimum  phase  velocities  of  the 
extensional  wave. 

We  begin  by  developing  the  equations  of  motion  of  an  infinitely  long,  layered  rod  with  a 
transversely  isotropic  core  and  a  transversely  isotropic  casing  using  the  three-dimensional 
equations  of  elasticity.  The  equations  of  motion,  which  are  general  in  axial  wavenumber  k  and 
circumferential  wavenumber  n,  are  then  specialized  to  the  case  of  the  n  =  0  axisymmetric  modes 
of  wave  propagation.  Simplifications  to  these  equations  for  the  case  of  material  isotropy  are 
presented. 

A  model  of  the  maximum  phase  velocity  of  the  extensional  waves  in  composite  isotropic 
rods  is  discussed.  It  is  found  that,  for  a  given  set  of  material  parameters,  this  maximum  phase 
velocity  is  a  function  of  the  ratio  of  the  outer  radius  of  the  casing  and  the  radius  of  the  core. 

The  extensional  wave  is  dispersive,  and  as  wavenumber  increases  to  infinity  the 
extensional  wave  phase  velocity  asymptotically  approaches  a  minimum.  This  minimum  velocity 
depends  on  the  material  parameters  of  the  core  and  casing  but  is  independent  of  system  geometry. 
The  ratio  of  the  maximum  phase  velocity  to  the  minimum  phase  velocity  is  the  range  of  dispersion 
of  the  system. 

Methodologies  for  controlling  either  the  maximum  phase  velocity  of  the  extensional  wave 
or  the  range  of  dispersion  in  composite  rods  of  isotropic  materials  are  presented.  The  frequency 
spectrum  of  the  extensional  wave  for  four  cases  of  an  isotropic  core  and  casing  are  developed,  and 
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the  dispersion  characteristics  are  discussed. 


Finally,  dispersion  of  a  composite  rod  consisting  of  an  isotropic  core  and  transversely 
isotropic  casing  is  investigated.  The  transversely  isotropic  casing  is  a  fiber-reinforced  composite 
material.  A  model  of  the  maximum  phase  velocity  of  the  extensional  waves  in  the  system  is 
developed  and  compared  to  the  model  for  composite  isotropic  rods.  Because  of  a  priori 
assumptions  of  the  dynamic  behavior  of  the  transversely  isotropic  casing,  the  asymptotic  minimum 
phase  velocity  of  the  extensional  wave  cannot  be  investigated. 

The  frequency  spectrum  of  the  extensional  wave  in  a  rod  with  an  isotropic  core  and 
transversely  isotropic  matrix  is  developed  and  compared  to  that  of  a  composite  isotropic  rod.  The 
isotropic  core  of  both  rods  is  identical.  The  transversely  isotropic  casing  consists  of  a  matrix 
material  that  is  identical  to  the  casing  of  the  composite  isotropic  rod.  The  geometries  of  each 
system  are  selected  so  that  the  maximum  phase  velocity  of  the  extensional  wave  in  each  rod  is 
identical.  The  characteristics  of  the  dispersion  of  the  extensional  wave  in  the  two  rods  are 
discussed. 

PREVIOUS  INVESTIGATIONS 

Beginning  with  the  work  of  Pochhammer  (1876)  and  Chree  (1889),  the  modes  of  wave 
propagation  in  a  homogeneous  isotropic  circular  rod  have  been  investigated  in  considerable  detail. 
A  concise  discussion  of  the  frequency  spectrum  of  circular  rods  can  be  found  in  Achenbach  (1987). 
Dispersion  of  axisymmetric  modes  in  composite  elastic  rods  has  been  investigated  by  McNiven  et 
al.  (1963),  Armenakas  (1965),  Mengi  and  McNiven  (1967),  McNiven  and  Mengi  (1967),  and  Lai 
(1971).  However,  in  all  of  these  investigations,  the  composite  rod  consisted  of  an  isotropic  core 
and  an  isotropic  cladding. 

A  limited  number  of  investigations  exists  in  the  area  of  wave  propagation  in  cylinders 
composed  of  anisotropic  materials,  particularly  based  on  the  three-dimensional  theory  of  elasticity. 
Chree  (1890)  may  have  been  the  first  investigator  of  axisymmetric  wave  propagation  in  an 
anisotropic  bar.  In  his  analysis,  the  equations  of  motion  of  a  transversely  isotropic  bar  of  any  cross 
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section  were  developed  and  then  specialized  for  the  case  of  a  circular  rod.  Approximate  solutions 
were  obtained  by  expanding  the  displacements  in  powers  of  the  radius  r.  Later,  Morse  (1954) 
developed  the  exact  solution  of  the  problem  studied  by  Chree  (1890)  and  showed  that  the  solution 
reduced  to  the  Pochhammer  solution  for  the  isotropic  case.  No  numerical  results  were  presented 
by  either  Chree  (1890)  or  Morse  (1954).  Mirsky  (1965a,  b)  investigated  axisymmetric  and 
nonaxisymmetric  wave  propagation  in  transversely  isotropic  cylinders  and  used  displacement 
potentials  to  solve  the  equations  of  motion.  All  the  above  investigations  were  for  a  rod  or  a  hollow 
cylinder,  but  not  a  layered  rod. 

To  the  author’s  knowledge,  the  case  of  wave  propagation  in  a  composite  rod  consisting  of 
an  isotropic  core  and  a  transversely  isotropic  casing  has  not  been  investigated  to  date. 

TRANSVERSELY  ISOTROIC  MATERIALS 

Many  natural  and  artificial  materials  are  transversely  isotropic.  For  example,  hexagonal 
crystals  (including  beryllium,  cadmium,  magnesium,  titanium,  and  zinc)  are  transversely  isotropic. 
Ice  is  also  transversely  isotropic. 

A  unidirectional,  fiber-reinforced  composite  material  is  transversely  isotropic.  This  type  of 
material  is  highly  anisotropic  because  the  stiffness  and  strength  in  the  fiber  direction  are  of  the 
order  of  the  values  for  the  fiber,  and  are  thus  very  high,  while  the  stiffness  and  strength  transverse 
to  the  fiber  direction  are  of  the  order  of  the  values  for  the  matrix,  and  are  thus  much  lower.  Carbon 
and  graphite,  which  are  often  used  in  fiber-reinforced  materials,  are  themselves  transversely 
isotropic. 

The  constitutive  behavior  of  transversely  isotropic  materials  is  characterized  by  five 
independent  elastic  moduli.  At  every  point  in  such  materials,  there  exists  a  plane  of  isotropy  in 
which  the  elastic  moduli  of  the  material  are  equal  in  all  directions.  The  unit  vector  normal  to  the 
special  plane  of  isotropy  in  a  transversely  isotropic  material  is  called  the  preferred  direction  of  the 
material.  Comprehensive  discussions  of  material  symmetry  and  elastic  moduli  can  be  found  in 
Jones  (1975),  Christensen  (1979),  and  Sokolnikoff  (1987). 
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WAVE  PROPAGATION  IN  COMPOSITE  BODIES  AND  COMPOSITE  MATERIALS 


A  composite  body  consists  of  two  or  more  materials  that  are  macroscopically  distinct  and 
are  joined  together  in  order  to  achieve  some  useful  property.  The  frequency  equation  of  a 
composite  body  is  developed  from  the  imposed  boundary  conditions.  If  there  are  no  a  priori 
assumptions  of  the  behavior  of  the  system,  which  restrict  investigations  (e.g.,  to  low 
wavenumbers),  the  frequency  spectra  of  the  modes  of  wave  propagation  in  a  composite  body  can 
be  developed  for  all  wavenumbers.  This  is  not  true  for  wave  propagation  in  a  composite  material. 

In  a  composite  material,  it  is  possible  to  define  a  representative  volume  element  (Hashin, 
1983).  Representative  volume  elements  are  large  compared  to  the  dimensions  of  the  components 
of  the  material.  A  necessary  characteristic  of  a  composite  material  is  statistical  homogeneity.  In  a 
statistically  homogeneous  composite  material,  all  global  characteristics  (such  as  volume  fractions, 
two-point  correlations,  etc.)  are  the  same  in  any  representative  volume  element,  irrespective  of  its 
position  within  the  material. 

The  ejfective  elastic  constants  of  a  composite  material  are  used  to  define  the  relations 
between  the  averages  of  field  variables,  such  as  stress  and  strain.  The  classical  approach  in  the 
analysis  of  homogeneous  composite  materials  is  to  assume  that  the  field  equations  of  elasticity  are 
valid,  with  the  effective  properties  of  the  composite  material  replacing  the  usual  elastic  properties 
of  a  homogeneous  continuum. 

Wave  propagation  in  deformable  composite  materials  can  be  broadly  divided  into  two 
categories  (Christensen,  1979).  If  the  wavelength  of  the  response  of  the  material  is  long  compared 
to  the  scale  of  the  representative  volume  element,  the  material  response  is  governed  by  the  effective 
properties  of  the  material.  In  this  case,  the  equations  of  motion  for  the  composite  material  are 
identical  to  those  for  honiogeneous  materials.  If  the  wavelength  of  the  response  is  not  ideally  long 
with  respect  to  the  representative  volume  element,  very  complicated  dynamic  effects  occur.  These 
effects  include  wave  reflection  and  refraction  at  the  interfaces  of  the  components  of  the  composite 
material.  These  phenomena  are  most  important  in  the  field  of  ultrasonics. 
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In  this  investigation,  we  assume  that  the  dynamic  behavior  of  the  composite  material  is 
governed  by  the  effective  properties.  As  a  result,  the  range  in  wavenumber  for  which  we  can 
develop  the  frequency  spectra  is  limited. 


5/6 

Reverse  Blank 


ANALYTICAL  FORMULATION 


SYSTEM  GEOMETRY  AND  APPROACH 

The  system  under  consideration  is  shown  in  figure  1.  It  consists  of  an  infinitely  long 
composite  rod  with  a  solid  core  of  radius  a  and  a  casing  of  outer  radius  b.  The  system  is  assumed 
to  be  linear  so  that  the  linearized,  three-dimensional  stress  equations  of  motion  can  be  used.  The 
system  displacements  and  stresses  are  defined  by  the  cylindrical  coordinates,  r,  0,  and  z.  The 
cylindrical  coordinate  directions  are  designated  as  the  1,  2,  and  3  directions,  respectively. 

We  begin  by  assuming  that  the  core  and  casing  are  composed  of  differing  transversely 
isotropic  materials.  Each  material  has  a  preferred  material  direction  collinear  with  the  longitudinal 
axis  of  the  cylinder.  The  equations  of  displacement  and  surface  stress  are  presented  as  functions  of 
geometry  and  material  parameters.  Simplifications  of  these  equations  (due  to  specializing  either 
the  core,  the  casing,  or  both,  to  the  case  of  isotropy)  are  discussed. 

The  frequency  equation  of  the  system  is  then  developed  from  the  assumptions  of 
(1)  continuity  of  displacements  and  stresses  at  the  interface  between  the  core  and  casing  and  (2)  a 
traction-free  boundary  condition  on  the  surface  of  the  rod. 


DISPLACEMENTS  AND  STRESSES  IN  A  COMPOSITE  ROD 

Solutions  to  the  equations  of  motion  of  transversely  isotropic  cylinders  were  developed  in 
a  previous  report  (Berliner,  1995).  The  resulting  displacements  and  stresses  are  repeated  here  but 
in  a  more  general  form  to  allow  for  systems  consisting  of  multiple  layers. 

The  displacement  components  are 
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where  the  subscript  j  indicates  the  components  and  properties  pertaining  to  layer y;  n  is  the 
circumferential  wavenumber;  i  =  ri  ;  k  is  the  axial  wavenumber,  co  is  the  radian  frequency; 
and  are  Bessel  functions  of  order  n  of  the  first  and  second  kind,  respectively;  and  •/„  (Ci^)  ^ 
denotes  8/^  ( /dr ,  etc.  In  this  investigation,  /  =  1  is  used  to  indicate  properties  and 

components  pertaining  to  the  solid  core,  and  7  =  2  is  used  to  indicate  properties  and  components 
pertaining  to  the  casing.  Also,  for  finiteness  of  displacements  and  stresses  at  r  =  0 ,  the 
coefficients  Bjj ,  ,  and  must  be  equal  to  zero. 


The  parameters  in  equations  (1)  through  (6)  are  derived  from  the  properties  of  the  material 
as  follows: 
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and  where  jP  is  the  density  of  the  material  and  the  coefficients  ,  jC^2  ’  elastic 

constants  of  the  material  of  layer  j.  These  elastic  constants  are  related  to  the  more  familiar 
engineering  constants  as 
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jBj  =  the  transverse  bulk  modulus, 
pi  =  the  longitudinal  shear  modulus, 
p-f=  the  transverse  shear  modulus, 
pi  =  the  longitudinal  Young’s  modulus, 

=  the  transverse  Young’s  modulus, 
jVl  -  the  longitudinal  Poisson’s  ratio, 
jVf  =  the  transverse  Poisson’s  ratio. 

“Transverse”  refers  to  properties  in  the  plane  of  isotropy  and  “longitudinal”  refers  to  properties 
along  the  material  preferred  direction. 


Effective  Engineering  Constants  for  a  Fiber-Reinforced  Matrix  Material 

When  the  transversely  isotropic  material  is  a  fiber-reinforced  matrix  material, 

etc.,  are  the  effective  elastic  constants  of  the  material.  These  effective  elastic  constants  are  defined 
from  effective  engineering  constants.  For  this  investigation,  we  use  the  definitions  for  the  effective 
engineering  constants  derived  by  Jones  (1975)  as  follows: 
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where  £  ,  ,G„,  v.  o.  and  y  are  the  Young’s  modulus,  shear  modulus,  Poisson’s  ratio, 

density,  and  volume  fraction  of  the  matrix  material  of  layer  j,  respectively,  and  where 

j^f>  jGp  pp  Pp  and^Vy  are  the  Young’s  modulus,  shear  modulus,  Poisson’s  ratio,  density,  and 

volume  fraction  of  the  fiber  of  layer  j,  respectively.  The  transverse  bulk  modulus  of  each  layer  is 
defined  by  Hashin  (1979)  as 
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jEj 


L 
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Specialized  Case  of  Isotropy 

We  can  specialize  the  constitutive  behavior  of  the  core  or  the  casing  from  transversely 
isotropic  to  isotropic  by  reducing  the  number  of  elastic  constants  from  five  to  two  using  the 
following  relationships: 
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These  elastic  constants  are  written  in  terms  of  the  engineering  constants  (also  known  as 
technical  constants)  of  the  material  as  follows: 


(l  +  .v)  (1-2  .V)’ 


(l  +  .v)(l-2  .V)’ 


2(l+^.v)’ 


where  jE  is  the  Young’s  modulus  and  is  the  Poisson’s  ratio  of  the  material.  By  normalizing 

and  jC^2  fu  ’  elastic  constants  simplify  to 


2(l-.v) 

>^11  “  (1-2  V)  ^^44’ 

'  (15 

Substitution  of  equations  (13)  and  (15)  into  equations  (7)  and  (8)  results  in  the  following 
simplifications: 
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BOUNDARY  CONDITIONS 

To  complete  the  analytical  formulation  of  the  frequency  equation  of  the  system  requires 
consideration  of  the  boundary  conditions.  At  the  solid-solid  interface,  r  =  a ,  the  continuity 
conditions  require  that  the  three  displacement  components  and  the  three  surface  stress  components 
for  the  core  and  casing  be  equal  (Fung,  1977).  The  traction-free  boundary  condition  on  the  outer 
surface  of  the  composite  rod  requires  that  the  three  surface  stress  components  are  zero  at  r  =  b. 

These  boundary  conditions  can  be  written  as 

t iV’  l^zr>  l“r’  r  =  a  ~  ^2^’  2'^zr’  2'^r0’  2^r>  2^z^  r  =  a 

and 

[2'^rr’ 2^’ 2'^r0]  r  =  fc  ~  P’ d]  r  =  fc  •  (1^) 

Equations  (17)  and  (18)  represent  a  system  of  nine  linear  algebraic  equations  in  unknown 
wave  amplitudes  that,  in  matrix  form,  are 

[Lp,^]  {c}  =  {0}  (p  =  1,2,...9;<7=1,2,...9),  (19) 

where  [L^  is  a  9-by-9  matrix  whose  components  are  found  from  equations  (1)  through  (6)  and 


15 


where  {  c}  is  a  9-by-l  column  vector  of  the  unknown  amplitude  coefficients  Aj  i,  A21,  A31,  A12, 
Bi2,  A22,  B22,  A32,  and  B32.  The  components  of  [L^  are  defined  in  appendix  A. 

The  solution  of  equation  (19)  is  nontrivial  when  the  determinant  of  the  coefficients  of  the 
wave  amplitudes  {c}  vanishes,  i.e.,  when. 


K,\  =  0-  <20) 

Equation  (20)  is  the  frequency  equation  of  the  system.  When  the  material  and  geometrical 
parameters  are  given,  equation  (20)  yields  an  implicit  transcendental  relationship  between 
wavenumber  and  frequency. 
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NUMERICAL  ANALYSIS 


FREQUENCY  EQUATION  FOR  AXISYMMETRIC  MODES 

We  begin  by  specializing  the  coupled  system  for  the  case  of  the  n  =  0  axisymmetric 
modes  of  wave  propagation.  When  n  =  0 ,  the  displacements  of  the  system  are  independent  of  0 . 
The  motions  separate  into  torsional  modes  consisting  of  jUq  (r,  0,  z,  t)  only  and  into  longitudinal 

modes  consisting  of  u  (r,  0,  z,  t)  and  u  (r,  0,  z,  t) .  In  this  investigation,  we  restrict  our  results 
to  longitudinal  modes  of  the  system. 

When  n  =  0 ,  many  of  the  components  of  the  matrix  [L^  become  equal  to  zero,  and  as 
a  result,  equation  (20)  separates  into  the  product  of  the  two  subdeterminants 

Det^Det2  =  0,  (21) 

where 

=  h,  q\  (P  =  1’  2,  4,  6,  7,  S-q  =  1,  2,  3, 4. 7,  8)  (22) 

and 

=  h,q\  ^p  =  3,  5, 9-,q  =  5,  6, 9) .  (23) 

The  components  of  the  determinants  are  defined  in  appendix  A,  with  n  =  0. 

Equation  (21)  is  satisfied  if  either  Det^  orDet2  is  equal  to  zero.  The  case  Det-^  =0  is  the 

frequency  equation  for  axisymmetric  longitudinal  modes  in  a  composite  rod.  The  case  Det2  =  0 

is  the  frequency  equation  for  axisymmetric  torsional  modes  in  a  composite  rod.  Since  only  the 
longitudinal  modes  are  of  interest,  we  will  only  seek  solutions  to 

Derj  =  0.  (24) 
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SOLUTIONS  OF  THE  FREQUENCY  EQUATION  FOR  LONGITUDINAL  MODES 

The  arguments  of  the  Bessel  functions  in  equation  (24)  can  be  real,  imaginary,  or  complex. 
Therefore,  equation  (24)  is  a  complex  algebraic  equation.  Solution  of  this  equation  requires  that 
the  real  and  imaginary  parts  vanish  simultaneously.  For  any  wavenumber,  there  are  an  infinite 
number  of  frequencies  that  are  solutions  to  equation  (24).  Typically,  root  searches  for  each  branch 
of  the  spectrum  begin  at  a  cutoff  frequency.  The  wavenumber  is  then  increased  by  a  small 
increment  and  the  search  for  the  next  root  is  conducted  within  the  neighborhood  of  the  previous 
root.  This  process  is  repeated  over  the  wavenumber  band  of  interest.  Spurious  roots  can  also  occur 
at  the  frequencies  where  arguments  of  the  Bessel  functions  vanish.  These  roots  can  be  determined 
in  advance  and  eliminated  from  the  frequency  spectrum.  The  root  searches  were  conducted  using 
Mueller’s  method  (Press  et  al.,  1986).  Mueller’s  method  is  an  algorithm  that  finds  a  real  or  complex 
zero  of  a  complex  function,  which  in  this  case  is  equation  (24). 
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COMPUTATIONAL  RESULTS 


The  extensional  wave  is  the  lowest  and  is  commonly  accepted  to  be  the  most  important 
longitudinal  mode  (Achenbach,  1987;  McNiven  and  Mengi,  1967).  This  mode  is  characterized  by 
primarily  axial  displacements.  In  this  investigation,  we  limit  the  following  computational  results 
to  the  extensional  wave. 

These  results  are  discussed  as  functions  of  normalized  frequency  Q ,  normalized 
wavenumber  6 ,  normalized  phase  velocity  c ,  and  the  ratio  of  b  to  a,  where  these  normalized 
variables  are  defined  as 


a  =  radius  of  core, 
b  =  outer  radius  of  casing. 


Q  =  — ,  normalized  frequency, 
l^B 

5  =  ka,  normalized  wavenumber, 

c  =  normalized  phase  velocity, 
o 


and  where 


(25) 


is  the  bar  wave  velocity  in  the  core  of  the  composite  rod  (Achenbach,  1987).  Using  the  root  search 
technique  described  above,  the  frequency  spectrum  was  developed  over  real  normalized 
wavenumbers  from  0  to  0.1. 

ISOTROPIC  CORE  AND  ISOTROPIC  CASING 

The  first  cases  to  be  investigated  are  composite  rods  consisting  of  an  isotropic  silica  glass 
core  and  linear  isotropic  elastomeric  casings.  (Note  that  the  material  properties  for  all  cases 
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investigated  can  be  found  in  appendix  B.) 


Figure  2  is  the  frequency  spectra  of  the  extensional  wave  in  a  composite  rod  with  a  rubber 
casing  and  five  values  of  R.  Two  comments  can  be  immediately  made  about  the  behavior  of  these 
systems  as  functions  of  R  and  of  5:  (1)  As  increases,  the  normalized  phase  velocity  c  decreases 
over  all  normalized  wavenumbers  and  (2)  For  any  R,  the  normalized  phase  velocity  as  5  ->  0  is 
different  from  the  normalized  phase  velocity  as  5  ->  <» .  Because  the  phase  velocity  of  the 
extensional  wave  is  not  constant  with  wavenumber,  it  is  dispersive. 

The  larger  the  range  of  phase  velocities,  the  more  dispersive  is  the  extensional  wave. 
McNiven  and  Mengi  (1967)  established  the  extent  of  this  range,  as  functions  of  geometry  and 
material  properties,  by  defining  the  velocities  bounding  the  frequency  spectrum  of  this  mode. 

Maximum  Phase  Velocity  of  the  Extensional  Wave 

Using  the  classical  kinematic  assumption  that  plane  cross  sections  remain  plane  during 
deformation,  McNiven  and  Mengi  (1967)  found  that  the  maximum  phase  velocity  associated  with 
the  extensional  wave  occurred  when  k  =  0  and  that  this  velocity  was  equal  to  the  bar  wave 
velocity  of  the  composite  rod.  They  defined  this  velocity  as 


where 


and 


E  A+E^\ 
1  1  2  y  j 

A 


E  A^^E^ 
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In  equation  (28),  Ej^  is  the  effective  Young’s  modulus  of  the  composite  rod,  is  the  effective 

density  of  the  rod,  is  the  cross-sectional  area  of  the  core,  2^  is  the  cross-sectional  area  of  the 

casing,  and  A  is  the  total  cross-sectional  area  of  the  rod.  Equations  (27)  and  (28)  are  sometimes 
called  a  “rule  of  mixtures”  (Jones,  1975).  Lai  (1971)  later  derived  equation  (26)  from  the  exact 
frequency  equation  by  using  low  argument  approximations  for  the  Bessel  functions. 

We  can  put  equation  (26)  in  a  more  convenient  form  by  rearranging  and  substituting 
2  [22) 

jA  =  Tta  ,  ^  =  K\b  -a  l,andJ?  =  (b/a),  resulting  in 


(29) 


Inspection  of  equation  (29)  shows  that  the  maximum  phase  velocity  of  the  extensional  wave  in  the 
composite  rod  is  bounded  by  the  phase  velocity  of  the  bar  wave  in  the  core  as  /?  — >  1  and  by  the 
phase  velocity  of  the  bar  wave  in  the  casing  as  R  — >00,  where  the  bar  wave  velocity  in  the  casing 
is  defined  by 


In  the  case  of  a  silica  glass  core  with  a  rubber  casing,  the  phase  velocity  of  the  bar  wave  in 
the  casing  is  less  than  the  phase  velocity  of  the  bar  wave  in  the  core.  Therefore,  ^^E/ is  the 

upper  bound  and  J^E/^  is  the  lower  bound  of  the  maximum  phase  velocity  of  the  extensional 
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wave  in  the  composite  rod,  and  the  phase  velocity  decreases  as  R  increases,  as  was  noted  in  figure 
2.  Figure  3  shows  the  same  frequency  spectra  as  figure  2  but  only  plotted  to  6  =  0.01  and 
Q.  =  0.01  to  emphasize  the  progressive  reduction  in  phase  velocity  with  increasing  and  the 
nearly  linear  behavior  of  the  extensional  wave  at  low  wavenumbers. 


Minimum  Phase  Velocity  of  the  Extensional  Wave 

The  composite  rods  in  this  section  represent  composite  bodies  composed  of  two  isotropic 
materials.  As  discussed  earlier,  since  the  frequency  equation  for  this  system  is  exact,  the  frequency 
spectra  of  the  extensional  wave  can  be  developed  over  all  wavenumbers.  Therefore,  we  can  analyze 
the  behavior  of  these  systems  as  ^  0  (as  was  done  above)  and  the  asymptotic  behavior  of  these 

systems  as  k—^^.  Mengi  and  McNiven  (1967)  found  that  all  axisymmetric  waves  in  an  infinitely 
long,  composite  rod,  consisting  of  an  isotropic  core  and  an  isotropic  casing,  asymptotically 
approach  the  phase  velocity  of  certain  classical  waves  as  ^  >  oo .  These  classical  waves  are 
Rayleigh  waves  in  the  casing,  Stonely  waves  at  the  interface  of  the  core  and  casing,  and  shear 
waves  in  the  core  or  casing. 

Composite  rods  of  this  geometry  can  be  grouped  into  one  of  four  classes.  First,  the  rods  are 
grouped  into  one  of  two  classes  determined  by  whether  the  shear  wave  velocity  in  the  casing  is 
faster  or  slower  than  the  shear  wave  velocity  in  the  core.  The  material  in  which  the  shear  wave 
velocity  is  less  is  called  the  acoustically  denser  of  the  two  materials.  These  two  classes  are  each 
further  subdivided  into  two  classes,  depending  upon  the  ability  of  the  composite  rod  to  propagate 
Stonely  waves  at  the  interface.  A  brief  description  of  these  four  classes,  as  well  as  the  asymptotic 
phase  velocity  of  the  extensional  wave  for  each  class,  is  repeated  in  table  1  to  aid  the  reader. 
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Figure  3.  Low  Wavenumber-Frequency  Spectra  of  Extensional  Waves  in  Composite  Rods  With 
a  Silica  Glass  Core  and  a  Rubber  Casing  as  a  Function  of  R 


Table  1.  Classes  of  Composite  Rods 


Class 

Description 

I 

This  rod  has  an  acoustically  dense  core  and  the  system  propagates 
Stonely  waves.  The  extensional  wave  velocity  asymptotically 
approaches  the  Stonely  velocity. 

II 

This  rod  also  has  an  acoustically  dense  core  but  does  not  propagate 
Stonely  waves.  The  extensional  wave  velocity  asymptotically 
approaches  the  shear  wave  velocity  in  the  core. 

III 

The  casing  of  this  rod  is  acoustically  denser  and  the  system 
propagates  Stonely  waves.  The  extensional  wave  asymptotically 
approaches  the  velocity  of  Rayleigh  waves  in  the  casing. 

IV 

This  rod  also  has  an  acoustically  dense  casing  but  does  not  propagate 
Stonely  waves.  However,  the  extensional  wave  velocity  in  this  type  of 
rod  also  asymptotically  approaches  the  velocity  of  Rayleigh  waves  in 
the  casing,  as  it  does  in  the  Class  III  type  of  rods. 

The  frequency  equation  for  Stonely  waves  is  determined  by  the  boundary  conditions  at  the 
interface  of  two  materials  (Achenbach,  1987).  Stonely  waves  will  propagate  at  the  interface  of  two 
materials  only  when  the  roots  of  this  frequency  equation  are  real.  When  Stonely  waves  do  exist, 
they  are  nondispersive. 


The  composite  rod  of  a  silica  glass  core  and  a  rubber  casing  is  either  Class  III  or  Class  fV 
because  the  shear  wave  velocity  in  the  casing  is  slower  than  the  shear  wave  velocity  in  the  core. 
As  noted  above,  the  extensional  wave  velocity  of  either  of  these  classes  of  composite  rods 
asymptotically  approaches  the  velocity  of  Rayleigh  waves  in  the  casing,  regardless  of  whether  the 
system  propagates  a  Stonely  wave  or  not.  Therefore,  we  do  not  have  to  investigate  the 
characteristics  of  the  roots  of  the  frequency  equation  for  Stonely  waves  in  this  system. 

The  Rayleigh  velocity  in  the  casing  is  approximately  given  by 

(0.862+  1.14  2V) 

where 
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is  the  shear  wave  velocity  in  the  casing.  Equation  (31)  represents  the  minimum  phase  velocity  of 
the  extensional  wave  in  the  composite  rod.  Contrary  to  the  maximum  phase  velocity  of  the 
extensional  wave,  the  minimum  phase  velocity  of  this  system  is  independent  of  R.  It  is  also 
interesting  to  note  how  this  asymptotic  behavior  differs  from  that  of  a  solid,  one-material  rod.  For 
that  case,  Achenbach  (1987)  has  shown  that  the  phase  velocity  asymptotically  approaches  from 
below  the  velocity  of  Rayleigh  waves  and  that  at  some  intermediate  wavenumber  the  phase 
velocity  has  a  minimum  value  less  than  . 


Controlled  Dispersion  Characteristics  of  the  Extensional  Wave 

McNiven  and  Mengi  (1967)  defined  a  measure  of  the  dispersion  of  the  extensional  wave  in 
a  composite  rod  as  the  ratio  of  the  two  bounding  velocities.  This  measure  of  dispersion  is  written  as 

2 

(33) 

where  Cq  is  the  maximum  phase  velocity  and  is  the  minimum  phase  velocity.  In  the  case  under 

investigation,  equation  (29)  defines  the  maximum  phase  velocity  and  equation  (31)  defines  the 
minimum  phase  velocity,  resulting  in 


Equations  (29),  (31),  and  (34)  can  be  used  to  control  the  dispersion  characteristics  of  the 
extensional  wave  in  composite  rods,  as  discussed  below. 
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Controlling  the  Maximum  Phase  Velocity  of  the  Extensional  Wave.  When  the  bar  wave 
velocity  of  the  casing  is  less  than  that  in  the  core,  the  maximum  phase  velocity  of  the  extensional 
wave  in  the  composite  rod  is  bounded  by  the  bar  wave  velocity  of  the  core  (upper  limit)  and  the 
bar  wave  velocity  of  the  casing  (lower  limit).  For  a  given  set  of  properties  for  the  core  and  casing, 
we  can  limit  the  maximum  phase  velocity,  within  these  bounding  velocities,  by  selecting  R. 
Rearranging  equation  (29),  we  have 


R  = 


fj 


(35) 


As  an  example  of  this  procedure,  we  arbitrarily  select  a  maximum  phase  velocity,  Cg ,  of 


1000  m/sec  and  determine  the  required  value  of  R  for  four  different  casing  materials,  designated 
as  “rubber,”  “lomod,”  “alcryn_low,”  and  “alcryn_high,”  each  with  silica  glass  as  the  core  material. 
As  mentioned  earlier,  the  properties  for  these  materials  are  listed  in  appendix  B.  When  the  material 
properties  and  =  1000  m/sec  are  substituted  into  equation  (35),  the  following  ratios  of  outer 

casing  radius  to  core  radius  are  found: 


rubber:  R  =  8.45, 
lomod:  R  =  9.37, 
alcryn_low:  R  =  8.79, 
alcryn_high:  R  =  8.46. 

Using  these  values  fori?,  along  with  the  pertinent  material  properties,  the  frequency 
spectrum  of  the  extensional  wave  was  developed  for  each  composite  rod.  Figure  4  is  a  plot  of  each 
of  these  spectra  along  with  the  line  representing  a  constant  normalized  phase  velocity  of  KKK)  m/ 
sec.  It  is  apparent  that,  although  the  maximum  phase  velocity  of  each  composite  rod  has  been 
limited  to  1000  m/sec,  the  range  of  dispersion  is  not  the  same.  This  is  because  the  minimum  phase 
velocity  of  each  composite  rod  is  different. 


We  use  equation  (31)  to  calculate  the  minimum  phase  velocity  of  each  composite  rod  with 
the  following  results: 
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Frequency  Spectra  of  Extensional  Waves  in  Composite  Rods  With 
Sihca  Glass  Core  and  Casings  of  Four  Different  Materials 


rubber:  2^^  =  68.19  m/sec, 
lomod:  2^^  =  251.61  m/sec, 
alcrynjow:  2^^^  =  167.98  m/sec, 
alcryn_high:  2^^  =  71.31  m/sec. 

The  range  of  dispersion  can  now  be  determined  by  substituting  Cq  =  Cg  =  1000  m/sec  and 
=  2Cr  for  of  the  materials  into  equation  (33)  with  the  following  results: 

rubber:  //  =  215.06, 
lomod:  //=  15.8, 
alcryn_low:  //=  35.44, 
alcryn_high:  H=  196.65. 

The  larger  the  range  of  dispersion,  H,  the  more  dispersive  is  the  composite  rod  (McNiven 
and  Mengi,  1967)  and  the  lower  in  wavenumber  is  the  phase  velocity  of  the  extensional  wave  when 
it  begins  to  deviate  from  the  maximum.  Inspection  of  the  above  calculations  for  H  indicate  that  the 
composite  rod  with  a  casing  of  lomod  material  is  the  least  dispersive  and  the  composite  rod  with  a 
casing  of  rubber  is  the  most  dispersive.  This  agrees  with  the  results  plotted  in  figure  4. 

Figure  5  shows  the  same  frequency  spectra  of  figure  4  but  plotted  over  a  smaller 
wavenumber  and  frequency  band.  This  figure  emphasizes  the  nearly  linear  behavior  of  the 
extensional  wave  at  the  low  wavenumbers  and  frequencies.  The  exact  phase  velocity  of  the 
extensional  wave  at  5  =  0.01  in  each  composite  rod  with  the  listed  material  as  the  casing  rhaterial 
is 


rubber:  c  =  717  m/sec, 
lomod:  c  s  973  m/sec, 
alcrynjow:  c  =  936  m/sec, 
alcryn_high :  c  =  7 34  m/sec . 

These  results  also  confirm  that  the  rod  with  the  casing  of  lomod  is  the  least  dispersive  and  that,  at 
6  =  0.01,  the  phase  velocity  of  the  extensional  wave  is  only  approximately  2.7  percent  less  than 
the  maximum  phase  velocity  of  1000  m/sec.  On  the  other  hand,  the  rod  with  the  casing  of  rubber 
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Figure  5.  Low  Wavenumber-Frequency  Spectra  of  Extensional  Waves  in  Composite  Rods 
With  a  Silica  Glass  Core  and  Casings  of  Four  Different  Materials 


is  the  most  dispersive  and,  at  6  =  0.01 ,  the  phase  velocity  of  the  extensional  wave  is  over  28 
percent  less  than  the  maximum  phase  velocity. 

Controlling  the  Range  of  Dispersion  of  the  Extensional  Wave.  Instead  of  controlling  the 
maximum  phase  velocity  of  the  extensional  wave  in  the  composite  rod,  we  may  wish  to  control  the 
range  of  dispersion.  We  can  use  equation  (34)  to  solve  for  R  for  any  desired  H  when  given  the 
material  properties  of  the  core  and  the  casing.  For  example,  if  the  casing  material  is  rubber,  and  if 
the  desired  values  of  H  are  10  or  100,  the  following  values  of  R  are  required: 

H  =  m-,R  =  47.8, 

//  =  100;  =  12.7. 

If  the  casing  material  is  lomod,  the  following  values  of  R  are  required  for  the  same  values  of  H: 

H  =  \0,R  =  12.9, 

H  =  m\R  =  3.42. 

Inspection  of  equation  (34)  shows  that  for  any  given  material  properties  of  the  core  and 
casing,  there  is  a  lower  bound  on  //  as  — > «» .  The  lower  bound  for  the  composite  rod  with  a 

rubber  casing  is  //  =  3.22,  and  the  lower  bound  for  the  composite  rod  with  a  casing  of  lomod  is 
H  =  3.16. 

ISOTROPIC  CORE  AND  TRANSVERSELY  ISOTROPIC  CASING 

The  next  case  to  be  investigated  is  a  composite  rod  consisting  of  an  isotropic  silica  glass 
core  and  a  transversely  isotropic  casing.  The  transversely  isotropic  casing  consists  of  a  fiber- 
reinforced  rubber  matrix  with  a  preferred  material  direction  collinear  with  the  longitudinal  axis  of 
the  rod.  The  matrix  material  has  the  same  material  properties  as  in  the  case  of  the  composite  rod 
with  an  isotropic  rubber  casing.  The  fibers  are  nylon  and  the  fiber  volume  fraction  is  10  percent. 
(The  material  properties  of  the  matrix  material  and  the  nylon  fibers  are  found  in  appendix  B.) 

When  the  material  properties  of  the  matrix  material  and  of  the  reinforcing  fibers,  along  with 
a  10-percent  volume  fraction  of  fiber,  are  substituted  into  equations  (11)  and  (12),  the  effective 
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engineering  constants  of  the  transversely  isotropic  casing  are  found  to  be 


2^7  =  1.54  X  lo’  N/m^, 

=  5.75  X  10®  N/m^ 

^Gj  =  5.75  X  10®  N/m^, 

2^^  =  5.64  X  10^  N/m^, 

=  1.67  X  10^  N/m^, 

2'^L  =  0.42, 

2V  j.  =  0.45, 

2p  =  1014kg/m^ 

Substitution  of  the  effective  engineering  constants  into  the  definitions  of  the  elastic 
constants  (see  equation  (10))  results  in 

Cl  1  =  2.12  X  lo”^  N/m^, 

C12  =  9.69  X  10®  N/m^, 

Ci3  =  1.3  X  10^  N/m^ 

C33  =  5.74  X  10*  N/m^, 

C44  =  5.74  X  10®  N/m^, 

C44  =  5.74  X  10®  N/m^ 

Maximum  Phase  Velocity  of  the  Extensional  Wave 

In  this  section,  we  develop  a  low-wavenumber  model  of  the  maximum  extensional  wave 
velocity  in  this  system  based  on  the  following  assumptions: 

1.  We  use  the  same  classical  kinematic  assumption  of  McNiven  and  Mengi  (1967).  That  is, 
we  assume  that  plane  cross  sections  remain  plane  during  deformation  and  that  the  extensional  wave 
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velocity  at  low  wavenumbers  is  defined  by  equation  (26)  as 


2.  Since,  at  low  wavenumbers,  the  extensional  wave  is  characterized  by  primarily 
longitudinal  displacements,  we  assume  that  the  longitudinal  Young’s  modulus  2^1  and  the  density 

2p  are  the  governing  material  properties  of  the  casing  for  this  mode.  Therefore,  we  define  the 
effective  Young’s  modulus  and  the  effective  density  of  the  composite  rod  as 


and 


Pe//  - 


(36) 


(37) 


respectively. 

When  equations  (36)  and  (37)  are  substituted  into  equation  (26),  along  with  =  na^, 

2A  =  Kib^-a^),  and/?  =  /?/ a,  the  low-wavenumber  model  of  the  phase  velocity  of  the 
extensional  wave  becomes 
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The  form  of  equation  (38)  is  identical  to  equation  (29)  when  ^  is  replaced  by  2^^  •  As  in  the  case 

of  a  composite  rod  consisting  of  an  isotropic  core  and  casing,  the  maximum  phase  velocity  of  the 
extensional  wave  in  this  system  is  bounded  by  the  phase  velocity  of  the  bar  wave  in  the  core  as 
^  ^  1  and  by  the  phase  velocity  of  the  bar  wave  in  the  casing  as  R  —>oo^  where  the  bar  wave 
velocity  in  the  casing  is  defined  as 


This  bounding  velocity  depends  on  the  volume  fractions  and  properties  of  the  matrix  material  and 
reinforcing  fibers  since  from  equation  (1 1)  we  note 


■Er  =  V.E+V  E  . 

J  L,  J  f  J  f  ]  m  ]  m' 


(40) 


Minimum  Phase  Velocity  of  the  Extensional  Wave 

The  composite  rod  in  this  section  represents  a  composite  body  composed  of  an  isotropic 
core  and  a  transversely  isotropic  casing.  Because  the  transversely  isotropic  casing  consists  of  a 
composite  material,  assumptions  regarding  its  dynamic  behavior  have  been  made  a  priori.  That  is, 
since  we  have  assumed  that  the  material  behavior  of  the  casing  can  be  represented  by  effective 
properties,  we  have  limited  our  investigation  to  wavelengths  that  are  long  compared  to  the 
representative  volume  element. 
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As  a  result,  we  can  discuss  the  behavior  as  ^  >  0  but  not  the  asymptotic  behavior  as 
This  prevents  us  from  developing  a  model  of  the  minimum  phase  velocity  of  the 
extensional  wave  of  this  system. 

Controlling  the  Maximum  Phase  Velocity  of  the  Extensional  Wave 

Since  we  cannot  define  a  minimum  phase  velocity  of  the  extensional  wave,  we  cannot 
define  a  measure  of  the  dispersion  of  this  mode.  However,  using  equation  (38),  we  can  control  the 
maximum  phase  velocity  of  the  extensional  wave. 

First,  as  an  example  of  the  effect  of  the  reinforcing  fibers  on  the  phase  velocity  of  the 
extensional  wave,  we  substitute  the  values  of  ^E,  ,  2^1  >  and  2P  into  equation  (38),  along  with 

R  =  8.45 ,  and  find  that  the  resulting  phase  velocity  is  1228  m/sec.  This  value  of  R  resulted  in 

a  phase  velocity  of  1(X)0  m/sec  in  the  composite  rod  with  an  isotropic  rubber  casing.  Therefore,  the 
addition  of  the  reinforcing  fibers,  without  changing  the  value  of  R,  has  increased  the  phase  velocity 
by  23  percent. 

Now,  by  rearranging  equation  (38)  to 


;  ~  2^L  ~  iP  2P  j 


(41) 


and  substituting  the  desired  value  of  =  KXX)  m/sec ,  we  find  that  R  must  equal  12.5.  This 

means  that,  when  10-percent  reinforcing  fibers  are  added  to  the  casing,  the  value  of  ;R  must  increase 
by  45  percent  to  maintain  a  maximum  phase  velocity  of  the  extensional  wave  of  1000  m/sec. 


The  frequency  spectrum  for  the  extensional  wave  in  the  composite  rod  with  a  silica  glass 
core  and  fiber-reinforced  rubber  casing  is  shown  in  figure  6.  This  system  has  an  R  =  12.5  so  that 
the  maximum  phase  velocity  is  1000  m/sec.  Also  shown  in  figure  6  is  the  frequency  spectrum  of 
the  composite  rod  consisting  of  the  silica  glass  core  and  isotropic  rubber  casing  and  a  line 
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_  -Rubber  casing 

_ Fiber-reinforced  casing 

_ 1 000  m/sec  constant  veiocity 


Figi^e  6.  Frequency  Spectra  of  Extensional  Waves  in  Composite  Rods  With 
Silica  Glass  Core  and  a  Casing  of  Isotropic  or  Fiber-Reinforced  Rubber 


representing  a  constant  normalized  phase  velocity  of  1000  m/sec.  Inspection  of  figure  6  shows  that 
the  extensional  wave  in  the  rod  with  the  transversely  isotropic  casing  has  a  distinct  change  in  phase 
velocity  at  approximately  5  =  0.035 .  The  model  of  the  extensional  wave  in  the  composite  rod 
assumes  only  axial  displacements  occur  and  uses  the  longitudinal  Young’s  modulus  to  define  the 
wave  velocity.  These  assumptions  are  valid  in  the  low  wavenumbers,  as  is  evidenced  by  the 
agreement  with  the  1000-m/sec  constant  phase  velocity  line.  However,  as  wavenumber  increases, 
radial  motions  couple  into  this  mode  and  the  simple  model  is  no  longer  valid. 

In  the  case  of  the  rod  with  a  transversely  isotropic  casing,  2^j  «  2^L’  might 

expect  that  as  radial  displacements  couple  into  the  longitudinal  mode,  the  phase  velocity  would 
decrease.  In  the  case  of  the  rod  with  an  isotropic  casing,  the  Young’s  modulus  is  constant  in  all 
directions.  Therefore,  we  might  expect  a  more  gradual  change  in  the  phase  velocity  as  the  radial 
displacements  begin  to  couple  into  this  mode.  However,  these  speculations  should  be  validated  by 
an  analysis  of  the  mode  shapes,  which  is  beyond  the  scope  of  this  investigation. 

Figure  7  shows  the  same  frequency  spectra  of  figure  6  but  plotted  over  a  smaller 
wavenumber  and  frequency  band.  This  figure  emphasizes  the  nearly  linear  behavior  of  the 
extensional  wave  at  the  low  wavenumbers  and  frequencies.  The  exact  phase  velocity  of  the 
extensional  wave  at  5  =  0.01  in  each  composite  rod  is 

isotropic  casing:  c  =  717  m/sec 
transversely  isotropic  casing:  c  =  875  m/sec. 

These  results  indicate  that  the  rod  with  the  transversely  isotropic  casing  is  less  dispersive  and  that 
at  5  =  0.01  the  phase  velocity  of  the  extensional  wave  is  approximately  12.5  percent  less  than  the 
maximum  phase  velocity  of  1000  m/sec.  The  rod  with  the  isotropic  rubber  casing  is  the  more 
dispersive,  and  at  5  =  0.01  the  phase  velocity  of  the  extensional  wave  is  over  28  percent  less  than 
the  maximum  phase  velocity. 


37 


Figure  7.  Low  Wavenumber-Frequency  Spectra  of  Extensional  Waves  in  Composite  Rods  With 
a  Silica  Glass  Core  and  a  Casing  of  Isotropic  or  Fiber-Reinforced  Rubber 


CONCLUSIONS 


The  maximum  phase  velocity  of  the  extensional  wave  in  composite  rods  with  isotropic 
cores  and  isotropic  casings  is  bounded  by  the  bar  wave  velocity  in  the  core  and  the  bar  wave 
velocity  in  the  casing.  With  proper  selection  of  the  ratio  of  the  outer  radius  of  the  casing  and  the 
radius  of  the  core,  the  maximum  phase  velocity  can  be  controlled  within  these  bounding  velocities. 
For  composite  rods  composed  of  isotropic  cores  and  isotropic  casings,  the  extensional  wave  phase 
velocity  has  an  asymptotic  minimum.  This  minimum  phase  velocity  can  be  used  to  calculate  the 
range  of  dispersion  of  the  extensional  wave.  For  the  types  of  systems  investigated  in  this  report, 
the  minimum  phase  velocity  is  the  Rayleigh  velocity  of  the  casing. 

A  low-wavenumber  model  of  the  maximum  phase  velocity  of  extensional  waves  in 
composite  rods  with  an  isotropic  core  and  a  transversely  isotropic  casing  was  developed.  It  was 
found  that  when  longitudinal  reinforcing  fibers  are  added  to  the  casing,  the  extensional  wave  phase 
velocity  increases.  The  maximum  phase  velocity  of  this  system  is  also  bounded  by  the  bar  wave 
velocity  in  the  core  and  the  bar  wave  velocity  in  the  casing,  where  the  bar  wave  velocity  of  the 
casing  is  defined  by  the  effective  longitudinal  Young’s  modulus  and  the  effective  density  of  the 
transversely  isotropic  material. 
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APPENDIX  A.  BOUNDARY  EQUATIONS 
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